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Abstract Tree-ring chronologies were developed from 
Abies mariesii Masters and Picea jezoensis ssp. hondoensis 
(Mayr) P. Schmidt collected at different altitudes of Mt. 
Norikura, central Japan. The tree-ring parameters of ring 
width and maximum density were measured by soft X-ray 
densitometry. The measurement series were detrended by 
fitting a 33-year cubic smoothing spline and autoregressive 
model. The correlation between species and between sites 
showed different responses of the species to climate in 
terms of ring width and maximum density. The correlation 
coefficient between sites within a species was higher than 
that between species for a site in the ring width, and the 
coefficient between species for a site was higher in maxi¬ 
mum density. The correlation coefficient between tree-ring 
chronology and monthly climate data set showed different 
responses of radial growth to climate. The different re¬ 
sponse was probably explained by the difference in the 
length of the growing season. High summer temperature 
increased the maximum density of the two species and the 
radial growth in Abies mariesii. Summer precipitation cor¬ 
related negatively to maximum density, but it did not affect 
the ring width of either species. The climatic signals could 
be extracted from ring-width and maximum-density chro¬ 
nologies of both species. 
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Introduction 

Tree growth is generally affected by climate changes, and 
many researchers have published papers on tree response to 
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climatic change and reconstruction of past climate using 
tree-ring chronology. 1-3 Many have focused on trees grow¬ 
ing under cold or dry climate 4,5 because the tree growth is 
expected to be more sensitive to climatic changes under a 
severe condition. In rather wet and warm region of temper¬ 
ate Asia, however, tree growth has been considered not to 
be as sensitive to climate factors. Thus, study sites were 
often chosen in high elevations or a northern area, expect¬ 
ing temperature to be a limiting factor of tree growth. 6-9 

The main islands of Japan are located in warm and cold 
temperate zones e.g. These geographical locations are not 
advantageous (for the above-mentioned reasons) to obtain 
good samples for dendroclimatological study, Recent 
work, however, is gradually extending the possibility of 
dendroclimatological study in Japan. Sweda 10 developed the 
ring-width chronology of Chameacyparis obtusa Endl. from 
Mt. Ontake, central Japan, and reconstructed past winter 
temperatures in the region. Yasue et al. 9 developed the 
chronologies of Picea glehnii Masters in terms of ring width 
and maximum density and discussed their responses to cli¬ 
mate in Hokkaido, northern Japan. D’Arrigo et al. 11 devel¬ 
oped ring-width chronology of Quercus dentata Thunb. in 
northeastern Hokkaido and indicated that the climatic sig¬ 
nal could be extracted from the chronology. Estimating the 
effect of climate on tree growth is essentially the basis for 
dendroclimatology. Only a few species have been surveyed 
in Japan for developing the chronology to discuss the effect 
of climate on tree growth. 

Natural forests of Abies spp. are widely distributed from 
the central to the northern subalpine region in Japan. 12 
Thus, the chronologies of these species present a good pos¬ 
sibility to establish a tree-ring network. Although Picea spp. 
are not distributed in as wide a range compared with Abies 
spp., they are the dominant species in the subalpine area. 
In this study, we developed chronologies of ring width and 
maximum density for Abies mariesii Masters and Picea 
jezoensis ssp. hondoensis (Mayr) P. Schmidt from different 
altitudes in central Japan. The purpose was to examine the 
potential of using these species for dendroclimatological 
study, including the response of radial tree growth to cli¬ 
mate changes. 
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Materials and methods 

Sites and sampled trees 

Mt. Norikura (3026 m above sea level) in central Japan was 
chosen for the sampling sites of this study (Fig. 1A). The 



Fig. 1. A Map showing the location of Mt. Norikura and Matsumoto 
City in which the meteorological station used in this study is located. 
B Sampling sites 


timberline is at 2450-2550m in altitude, and Japanese stone 
pine [Pinus pumila (Parras) Regel] is distributed above the 
timberline. Abies mariesii and Abies veitchii Lindl. are the 
most common tree species in subalpine natural forest 
around central Japan, and the former appears at higher 
altitudes than the latter. 13 In some cases, P. jezoensis ssp. 
hondoensis is mixed into the subalpine natural forests. 

Two sampling sites in mixed stands of A. mariesii and 
P. jezoensis ssp. hondoensis were selected at 2000 m 
(NRK4) and 2200 m (NRK5) in altitude on the east slope of 
the main ridge (Fig. IB). Because these two sites are located 
in natural forest, it is believed that no artificial disturbance 
has affected tree growth. The distance between the sites is 
so short (ca. 950 m on the map) there would be little differ¬ 
ence in climate throughout the year except for the tempera¬ 
ture difference due to altitude. In 1995 and 1997 two 
increment cores 5 mm in diameter were taken from each 
tree at breast height using an increment borer. Table 1 gives 
the details of the sites and sampled trees. 

Sample preparation and measurement 

The increment cores were mounted in the groove of wood 
sticks with vinyl acetate emulsion and cut into 1.56 or 
2.00 mm thick cross sections for soft x-ray densitometry. X- 
ray negatives were taken for air-dried sections after season¬ 
ing in a conditioning room (20°C, 60% RH). The nega¬ 
tive films were scanned using Dendro2003 (Walesch 
Electronic), and ring width and maximum density were 
measured. 

Crossdating 

Crossdating was performed by the following three steps. 
The first step was visual comparison of the densitometric 
pattern and synchronization of narrow rings (pointer years) 
among samples. The density fluctuation of the transition 
zone from earlywood to latewood gave effective informa¬ 
tion for crossdating. The second step was calculating the 
moving correlation. Correlation coefficients between each 
ring-width series of 50-year segments and the total mean 
series were calculated for every lag year from -80 to +80, 
and the lag years were listed when significant coefficients 
were obtained. Before calculation, the individual series of 
measurements were detrended by 33 years cubic spline 
filter 14 and prewhitened by fitting an autoregressive model. 
Steps 1 and 2 were attempted repeatedly until no lag in 
dating was detected visually or statistically. Finally, the 


Table 1. Summary of sites and sampled trees 


Code 

Site 

Location 

Altitude (m) 

Species 

No. of trees 

N4A 

NRK4 

36°06'58"N, 137°35'28"E 

1990-2000 

Abies 

5 

N4P 




Picea 

7 

N5A 

NRK5 

36°07'18"N, 137°34'59"E 

2180-2200 

Abies 

10 

N5P 




Picea 

14 


Abies, Abies mariesii Masters; Picea , Picea jezoensis ssp. hondoensis (Mayr) P. Schmidt. 
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crossdating was cross-checked by the program 
COFECHA, 15 which is widely used for statistical confirma¬ 
tion in dendrochronology and dendroclimatology. 

Chronology development 

To extract climatic signals, measurement series must be 
detrended to reduce variations due to age trends and natu¬ 
ral disturbances. 16 The tree-ring data were detrended by 
fitting the cubic smoothing spline with 33 years filter length 
by the program ARSTAN. 14 This filter reduces the ampli¬ 
tude to less than 50% in the range of 33 years and longer 
frequency. 17 Thus, the filter extracts the short-term varia¬ 
tion to examine the climatic signal. Because one of the 
purposes of this study was to confirm the possibility of 
extracting climatic signals from tree-ring chronologies, we 
believed that 33 years was the best filter length to reduce the 
influence of natural disturbance and age trend. The mean 
chronology was calculated by biweight robust means for 
removing the effect of extreme measurements on the mean 
value. 18 The standard chronologies were prewhitened by 
fitting an autoregressive model selected on the basis of the 
first minimum Akaike Information Criterion (AIC) to pro¬ 
duce residual chronologies. The residual chronologies are 
regarded as the chronologies in which the influence of past 
growth is removed. Thus, climatic signals are expected to be 
extracted more, as even long-term variation is lost from the 
chronologies. Both standard and residual chronologies con¬ 
sist of at least three cores for each year. 

Climate data 

The climate data sets (monthly mean temperature and 
monthly total precipitation, 1900-1990) observed at Matsu- 
moto Meteorological Station (36°15'06"N, 137°58'04"E, 
610 m in altitude) were selected for calculating simple corre¬ 
lation coefficients with tree-ring chronologies. The city of 


Matsumoto is located about 35 km eastward of the sampling 
sites (Fig. 1A). The monthly data set for the 14 months 
from September of the previous year to October of the 
current year from 1900 to 1990 was used to calculate the 
coefficients. 


Results and discussion 

Chronology 

Chronologies of ring width and maximum density were de¬ 
veloped for two species at two sites. Because the results of 
crossdating indicated that no absent or false ring exists in 
any of the measured cores, the cores were adopted for 
developing the chronologies. The length of chronologies of 
N4A, N5A, N4P, and N5P were 128 years (1867-1994), 195 
years (1802-1996), 173 years (1824-1996), and 159 years 
(1838-1996), respectively. The basic statistics for the stan¬ 
dard and residual chronologies are shown in Table 2. The 
mean sensitivity (MS), signal-to-noise ratio (SNR), and 
mean correlation between trees (MCT) were higher in re¬ 
sidual chronologies than in standard ones, which indicated 
that the common climatic signal was extracted more by 
fitting an autoregressive model because MS presents the 
variation due to yearly difference, and SNR and MCT are 
expressions of the strength of the observed common signal 
among trees in the ensemble. 17 Therefore, the residual chro¬ 
nologies were used for further analysis in this study (Figs. 2, 
3). The similarity of variation in four maximum-density 
chronologies seems to be high, as the increase/decrease of 
indices occur synchronously (e.g., in 1866,1922, and 1994). 
In contrast, the pointer years of ring width tend to appear 
differently between species (e.g., 1918, 1936, and 1990 in 
Abies and 1917, 1935, and 1991 in Picea). This result sug¬ 
gests that the response of radial growth is different between 
species even growing in the same sites. 


Table 2. Basic statistics of standard and residual chronologies 


Parameter 

Standard chronology 



Residual chronology 



MS 

SNR 

MCT 

FAC 

MS 

SNR 

MCT 

N4A 

Ring width 

0.11 

1.61 

0.24 

0.15 

0.12 

1.84 

0.27 

Max. density 

0.04 

0.68 

0.12 

0.19 

0.04 

0.67 

0.12 

N5A 

Ring width 

0.13 

3.67 

0.27 

0.35 

0.15 

4.12 

0.29 

Max. density 

0.05 

3.38 

0.27 

0.13 

0.05 

3.92 

0.30 

N4P 

Ring width 

0.12 

4.07 

0.37 

0.17 

0.15 

4.50 

0,39 

Max. density 

0.05 

2.27 

0.25 

0.16 

0.05 

2.55 

0.27 

N5P 

Ring width 

0.14 

9.85 

0.41 

0.40 

0.15 

10.61 

0.43 

Max. density 

0.05 

8.95 

0.39 

0.17 

0.05 

9.10 

0.39 


MS, mean sensitivity; SNR, signal-to-noise ratio; MCT, mean correlation between trees; FAC, 
first-order autocorrelation. 
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Calender year 


Fig. 2. Residual chronologies of ring width and sample depth. Arrows 
indicate pointer years. The minimum sample depth for each chronolo¬ 
gies is three cores 
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Calender year 

Fig. 3. Residual chronologies of maximum density. Arrows indicate 
pointer years 


Correlation between sites and between species 

The correlation coefficients were calculated between sites 
and between species for ring width and maximum density 
(Table 3). The length of the chronologies in the calculation 
was 91 years (1900-1990), in accord with the shortest chro¬ 
nology length of N4A (10 cores for the years 1896-1995). 

The correlation of ring width between two sites was 
highly significant within each species: 0.68 and 0.75 in N4A- 
N5A and N4P-N5P, respectively (Table 3). The coefficients 
between species within a site were somewhat smaller than 
between sites within a species. This result indicated that the 
growth response to environment was different in each spe¬ 


Table 3. Correlation coefficients between residual chronology from 
1900 to 1990 (91 years) 


Code 

N4A 

N5A 

N4P 

N5P 

N4A 


0.53 

0.60 

0.41 

N5A 

0.68 


0.70 

0.78 

N4P 

0.50 

0.54 


0.72 

N5P 

0.42 

0.53 

0.75 



All coefficients are significant at P < 0.01. Left bottom half shows 
correlation coefficients between the ring-width chronologies; right 
top half shows the correlation coefficients between maximum-density 
chronologies. 



Precipitation 






;. in, 
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Fig. 4. Correlation coefficients between chronologies and monthly cli¬ 
mate data. Dotted lines indicate the limit of significance (P < 0.05) 


cies, although the temperature conditions of their habitats 
estimated from the warmth index are similar. 19 Graumlich 20 
also reported that site-to-site differences are small com¬ 
pared with species-to-species differences in mixed conifer 
and deciduous forests in the region of relatively small cli¬ 
matic contrast. On the other hand, the coefficient calculated 
on maximum density between sites within a species was 
slightly lower than between species within a site. Thus, the 
factors affecting ring width and maximum density were 
probably different. 


Effect of monthly temperature on radial growth 

The correlation coefficients between tree-ring series and 
monthly climate data are shown in Fig. 4. The correlation 
coefficients between ring width and temperature during the 
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current growing season (July and August) were significant 
in two series of N4A and N5A. The positive correlation 
coefficients between ring width and temperature from April 
to August were also observed in N4P and N5P, but they 
were not significant. These results show that the response of 
radial growth to the climate is different between the species, 
and A. mariesii is more sensitive to summer temperature 
than P. jezonensis ssp. hondoensis in this area. Peterson and 
Peterson 21 reported that the relation between growth and 
current summer temperature varied in three species of 
subalpine fir ( Abies lasiocarpa (Hook. Nutt.), Engelmann 
spruce ( Picea engelmannii Parry), and subalpine larch 
(Larix lyallii Pari); growth of subalpine fir is positively 
correlated with July-August temperature, whereas that of 
Engelmann spruce is correlated with June-August tem¬ 
perature. Their results indicate that the duration affecting 
tree growth is longer in Picea than in Abies , an observation 
consistent with our results for P. jezoensis ssp. hondoensis 
and A. mariesii. 

There are some reports on the growing period of Abies 
spp in Japan. Kimura et al. 22 reported that Abies trees in Mt. 
Shimagare sprout their new needles in late June. Kaji 13 also 
observed that bud burst of A. mariesii occurred in the end of 
June in the subalpine area of Mt. Tengudake. Matsumoto 
and Negishi 23 reported that the photosynthetic capacity of 
A. veitchii increased immediately as an expansion of leaves 
in July and reached its maximum in August. Kajimoto (per¬ 
sonal communication) reported that the growth period for 
A. mariesii in northern Honshu is from the middle of June 
to August. These researchers indicated that the growing 
period of A. mariesii and A. veitchii in subalpine regions is 
July and August. 

The temperature during the same months positively 
affected the ring width of A. mariesii in our study. The 
significant correlation between ring-width chronologies and 
monthly temperature in July and August probably reflects 
the physiological activity in radial growth. The growing pe¬ 
riod for Picea spp. has been less well investigated compared 
with that for Abies. Chiba 24 reported that bud burst of P. 
jezoensis ssp. hondoensis occurs 2 weeks earlier than that of 
A. mariesii, planted at the same place in Hokkaido. Yasue 
et al. 7 found that initiation of cambial activity for Picea 
glehnii was at the beginning of May in northern Hokkaido. 
These studies suggested that the growing duration is likely 
to be longer for Picea spp. than for Abies spp. in the subal¬ 
pine regions of Japan. Therefore, the different growth re¬ 
sponse of the species observed in this study is believed to 
reflect species-specific growing characteristics, such as tim¬ 
ing of cambial reactivation. 

Effect of monthly temperature on maximum density 

The correlation between maximum density and monthly 
temperature was higher than that of ring width. The corre¬ 
lation coefficient is positively significant during early spring 
and the whole summer season except early summer. It was 
reported that maximum density correlated positively with 
summer temperature for Picea glehnii growing in northern 


Japan. 9,23 Because summer is the time for latewood 
formation, the positive influence of temperature is under¬ 
standable, as activated photosynthesis provides more pho- 
tosynthates to cambium and a differentiating zone. Warm 
spring probably stimulates physiological activity to burst 
buds earlier and induces a longer growing period. Conse¬ 
quently, a larger amount of carbohydrates is allocated for 
cell wall thickening throughout a growing season. Conkey 26 
discussed the contribution of spring conditions to the maxi¬ 
mum density of red spruce in Maine (United States). She 
assumed that early warming in spring permits early physi¬ 
ological activity, thereby providing increased photosyn- 
thates for a longer period of time during the growing 
season. Thus warm or early springs lead to high maximum 
density. 

A species-to-species difference is not clear for maximum 
density. Because maximum density is generally highly sensi¬ 
tive to summer climate in both Abies spp. and Picea spp., 27 
our result is interpreted as a response similar for the two 
species under similar microclimates at the sites. Therefore, 
it is thought that a site-to-site difference is larger than 
species-to-species difference for the maximum-density 
chronology. 

Effect of monthly precipitation on ring width and 
maximum density 

Precipitation significantly affected the maximum density 
negatively during the growing season for all sites and spe¬ 
cies and for ring width in December for Abies (Fig .4). 
Comparing the number of significant coefficients shows that 
precipitation had a much weaker influence on tree growth 
in these two sites. This fact suggests that moisture is not a 
limiting factor for tree growth in this area, but summer rain 
could reduce sunlight. The total precipitation at Matsumoto 
City is only about 1000 mm; the precipitation on the high 
mountain is more than 1000mm for the 3 months from July 
to September, which is measured at AMeDAS station 
(36°07'12"N, 137°33'30"E, 2730m) located at Mt. Norikura. 
Because the station works from mid-June to mid-October, 
the annual precipitation data are not available for estimat¬ 
ing the growth response over a year. Even though the total 
amount is unknown, the precipitation is expected to be 
enough for tree growth at the study sites. Sufficient precipi¬ 
tation through the growing season probably resulted in a 
low correlation between radial growth and precipitation, 
as the change of precipitation gives less impact on radial 
growth. Yasue et al. 9 reported that precipitation contributes 
rarely to the ring width and is not a limiting factor for Picea 
glehnii. 


Conclusions 

We developed the ring-width and maximum-density chro¬ 
nologies for Abies mariesii and Picea jezoensis ssp. 
hondoensis from two sites on Mt. Norikura. The species-to- 
species difference in ring width is larger than the site-to-site 
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difference in this area, and the species-specific response 
of radial growth to climate is inferred. On the other hand, 
maximum density is strongly affected by the micro¬ 
environment rather than the species characteristics. The 
ring width and maximum density of the two species are 
sensitive enough to temperature during the growing season 
that they are good tools for extracting common climatic 
signals. Thus, we conclude that both species used in this 
study are useful for further dendroclimatological study, and 
the development of longer chronologies promises successful 
results in the detection of long-term climate changes. 
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